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ABSTRACT 

We develop an optimized technique to extract density-density and velocity-velocity spectra 
out of observed spectra in redshift space. The measured spectra of the distribution of ha- 
los from redshift distorted mock map are binned into 2-dimensional coordinates in Fourier 
space so as to be decomposed into both spectra using angular projection dependence. With 
the threshold limit introduced to minimize nonlinear suppression, the decomposed velocity- 
velocity spectra are reasonably well measured up to scale k - 0.07 li Mpc"', and the measured 
variances using our method are consistent with errors predicted from a Fisher matrix analysis. 
The detectability is extendable to k ~ 0.1 li Mpc"' with more conservative bounds at the cost 
of weakened constraint. 
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1 INTRODUCTION 

The evolution of large scale structure, as revealed in the cluster- 
ing of galaxies observed in wide-deep redshift surveys has been 
one of key cosmological probes. Structure formation is driven by a 
competition between gravitational attraction and the expansion of 
space-time, which enables us to test our model of gravity at c osmo- 



logica l scales and the expansion of history of the Universe i Wan 
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Maps of galaxies where distances have been measured from 
redshifts show anisotr o pic deviations f rom the true galaxy dis- 
tribut i on (lYorkl l200d ; jpeacockl 1200 ll: ICollesj | 2003| - [Hawkins! 
20031 : IPercivall |2004 IZehavil l200a lUe Fevre 2005; TegmarM 
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recession velocities include 



Guzzo et al.l |2008[) . because galaxy 



components from both the Hubble flow and peculiar velocities. 
In linear theory, a distant observer should expect a multiplicative 
enhancement of the overdensity fi eld of tracers due to the pecu- 
liar motion along t he line of si ght i Davis & Peebles,. 1982; Kaisej 
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the observed redshift space distortion maps. However, it is not yet 
fully understood what the optimal technique is to practically de- 
compose the spectra as theory predicts. We propose a statistical 
technique to extract it up to the limit of theoretical estimation. 
Our method utilizes the distinct angular dependence of density- 
density and velocity-veclocity spectra to decompose them from 
two-dimensional redshift power spectra, and is consistent with the 
theoretical estimate from Fisher matrix analysis. 

We present the detailed formalism in the next section. The 
Fisher matrix analysis to decompose spectra is briefly reviewed, 
then we present the method to decompose spectra in an optimal 
way with mock data. We discuss statistical method to minimize the 
effect by nonlinear suppression. 



2 PECULIAR VELOCITY POWER SPECTRA 
EXTRACTION 

2.1 Theoretical Expectation of Decomposition Accuracy 

The observed power spectrum in redshift space is decomposed into 
spectra of density fluctuations and peculiar velocity fields in real 
space. The observed power spectra in redshift space, P, is given by. 



observed spectra in redshift space can be decomposed into both 
density-density a nd velocity-v elocity spectra using angular projec- 
tion dependence dSong & Percival 2009; Percival & White 200^ ; 
IWhite et al.ll2009l : ISong et alj|201oh . With a local linear bias, the 
real-space galaxy density field is affected, while the peculiar veloc- 
ity term is not. In this paper, we attempt to extract velocity-velocity 
spectra as an unbiased tool to trace the history of structure forma- 
tion. 

A theoretical formalism jWhite et al.ll2009l) was derived for 
forecasting errors when extracting velocity-velocity spectra out of 



P(k, fi, z) = {Ps,(k, z) + 2n-r(k) [Pjk, z)Pee(k, z)] 
+ fi*Pe,e,(k,z)\G(k,fi,cr,\ 



1/2 



(1) 



where Pgg is the galaxy-galaxy density spectrum, P@0 is the 
velocity-velocity spectrum (0 is the divergence of velocity map 
in unit of aH), and yu denotes the cosine of the angle between ori- 
entation of the wave vector and the line of sight. Because this de- 
composition is valid only at large scale and when the rotation of the 
velocity field is negligible, we focus on modes of < 0. 1 /jMpc"' 
jPueblas & Scoccimarrdl2009l) . The cross-correlation coefficient 
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Figure 1. Power spectta from mock map in 2D caitesian coordinate (k^^,kn). 



r{k) is defined as r(k) = Pge/ -^fPggPee- The density and velocity 
divergence are highly correlated for k < 0. 1 /iMpc~' so we a ssume 
that both are perfectly correlated, r(k) ~ 1 jWhite et alj2009l) . Then 
the density-velocity cross-spectrum becomes the geometric mean 
of the two a uto-spectra and we h ave only two free functions, Pgg 
and Pee. As IScoccimarrd j2004) clearly pointed out, the redshift 
space power spectrum is suppressed along line-of-sight due to the 
velocity dispersion of large-scale flow, and we follow his model by 
introducing a function G = exp(-^^/j^(Tj) where cr,, will be calcu- 
lated from linear theory. Considering the possibility that nonlinear 
dynamics, like Finger-of-Gods effect, might contaminate the power 
spectrum, we use this term to find a cut-off scale of fi to exclude 
data which cou ld be affected strongly by nonlinear dynamics. In- 
deed, l^mya^al. (2009J pointed out that tr,, calculated by linear 
theory does not match with result from N-body simulations if one 
tries to model the power spectrum at > 0.1 hMpc^' . This cut-off 
edge is defined by yu^ut = fth/^fr. where the value of a^, will 
be discussed later. 

We estimate the accuracy of decomposition of Pgg and Pqq 
from P using Fisher matrix analysis determining the sensitivity of a 
particular measurement. Fisher matrix for this decomposition, F'^'j^, 
is written as. 



^ dp, dpi, P(k,pf2(2nY 
where p„ = (Pgg, Pee)- The effective volume Vcff(^) is given by, 
nP 



(2) 



V,ff(P) ■■ 



nP+l 



(3) 



where n denotes galaxy number density. 

Derivative terms in Eq. (O are given by. 
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dPggik„Zj) 
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Pee(k„Zj) 
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Figure 2. The obsei-ved power spectra at scales, i=0.03, 0.05, 0.07 and 
O.OQ/iMpc"' (from top to bottom) are plotted with error bars at various fA. 
Solid curves are P(a{k,)i) (Kaiser effect alone) and dash curves are Pg, (/:,//) 
(including dispersion effect) from best fitting bias b(k). 



The diagonal elements of the inverse Fisher matrix indicate the esti- 
mated errors of decomposion accuracy. The variances of Pgg(kj,Zj) 
and P@@(ki,Zj) is given by. 



ar[Pgg(k„Zj)] 

o-[Peeik„z,)] 



^nr'(ki,zj) 

^|F'^%-'(k„ZJ). 



(5) 



2.2 2D power spectra from mock map 

We use the halo catalog ue from the time-stre aming mock map of 
the Horizon simulation iTevssier et al]|2008l) . and cut 1 (Gpc//i) 
cubic box at the median redshift z = 0.83, which contains 2.2 mil- 
lion halos. The fiducial cosmological parameters of the simulation 
are given by (fi,,, = 0.24, fi^ = 0,h = 0.7 2, as = 0.78, ng = 0.96 ) 
and the initial transfer function is given bv lSisenstein & Hul ( Il998h . 

The distribution of halos is modified according to their pecu- 
liar velocity to incorporate the redshift distortion effect. We adopt 
the distant observer approximation and measure the power spec- 
trum in (k^,kn) space. The density fluctuation field is constructed 
by assigning the halos to 512' grids for the fast Fourier transforma- 
tion (FFT) using the nearest grid point (NOP) method. Fig.[T]shows 
the resulting power spectrum. While linearly spaced bins in ik±,kn) 
are used in this plot for presentation purpose, we use bins in k and 
p for the following analysis, k is divided in Ak = 0.02 /iMpc"' lin- 
early equally spaced bins from k = 0.02 h Mpc"' to 0.2 h Mpc"' and 
/J is in 5 linear-bins from to 1 with equal spacing. The measured 
2D power spectra in (k,p) coordinate are shown in Fig.|2] 

The Gaussian variance is used to derive errors for each bin 
shown as error bars in Fig. |2] cr[P„i,(k, p)] = P(k, p) ^2IN{k, p) 
where N(k,p) is number of modes in Fourier space. We test this 
using an alternative method, jack-knife errors (we do not attemp to 
generate more samples as we are interested in mocking real observ- 
ables in a single patch). A total 64 jack-knife samples are prepared 
out of a single mock map by dividing each coordinate into 4 pieces. 
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Table 1. Best fitting biases b{k) at given scales k from k = 0.03 to 
0.09/iMpc"'. 



Both errors agrees well, and different bins weakly correlate with 
each other. 

Halo distribution is a biased tracer of the dark matter distribu- 
tion. Theoretical P,h(A:,/j) from Kaiser effect only is given by, 



(6) 



where Pmm(k) is the dark matter density-density spectra and h = 
b(k) is the halo bias for each given scale. Spectra P„,„,(k) and P@@(k) 
are given from the cosmological parameters used for the simula- 
tion, and the halo cross-correlation parameter rj, is set to be unity. 
It has been tested that r for dark matter-© is perfectly correlated at 
linear scales ^ < 0. 1 /7Mpc"' from simulation. Unfortunately, the 
same sanity check is not applicable for halo maps due to the insuffi- 
cient number of halo in each grid for direct velocity power spectra. 
Instead, the theoretical Pth(k,fi) is derived based upon rh(fe) = 1, 
and the possible departure from the unity is dectectable from mea- 
sured Pei&(k) at linear scales. 

The tracer bias is assumed not to be determined by theoret- 
ical formalism or by other experiment. Instead of applying scale 
independent bias, b(k) is varied independently for each ^-bin. We 
fit b(k) for each mode to get Pth(k, /j) (solid curves in Fig. O. 
In Table □ the best fit b{k) is given with 1-cr confidence level. 
Theoretical with fitted b(k) is over-plotted with the mea- 

sured Pob(k,fi) from the simulation in Fig. |2] We cut out scales 
k < 0.03 /iMpc ' due to our limited box size and k > OA /iMpc"' 
due to non-linear effects. 

Using Pth{k,fi), theoretical errors are estimated from Fisher 
matrix analysis. Un-fiUed black contours in Fig. [3] represent the 
theoretical expectation around b{kYP„,„,(k) and Po^ik). As it is 
prediced from halo bias model, measured bias is nearly scale in- 
dependent. 



2.3 Practical approach to extract peculiar velocity spectra 

Spectra Pgg{ki) and P@@(ki) are fitted simultaneously to Pa\,{ki,^p) 
where i and p denote k and /j bins respectively. Bias is not pa- 
rameterized to fit P„\y(ki,^p), instead, we use Pg^(kj). The fitting 
Pfn{ki,Hp) is given by 



Pu(ki,Hp) 



[Pggikd + 2^1 ^jPggPe,e, + ^i^Peeikdl 
G(k,iJj,cr„). 



(7) 



We consider the velocity dispersion effect from one-dimensional 
velocity dispersion cr,, which is given by. 



(8) 



This formula needs P^q which is what we want to measure. We 
will discuss how we calculate this term in the next paragraph. Eq.|7] 
is expected to be invalidated beyond some threshold. The observed 
modes are cut out when it goes beyond given the threshold limit cr,!, 
as kj/icutCv > CTth- The fiducial value is cr,i, = 0.24 which represents 
confidence of theoretical prediction up to 6% drop of G{kj,ficut, cry) 
from unity. 

The most important factor in the integration Eq.[8]is the ampli- 
tude of ^00, as scale-dependent factor of is tightly constrained 
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Figure 3. Contour plots are shown for decomposed Pgg(ki) and P00(kj) at 
k = 0.03, 0.05, 0.07 and 0.09 /iMpc"'. Unfilled black contours represent 
theoretical prediction from Fisher matrix analysis, and filled blue contours 
represent measured Pgg(ki) and PooC^:,) from mock map in redshift space. 



by CMB physics. The shape of the power spectra is determined be- 
fore the epoch of matter-radiation equality. When the initial fluc- 
tuations reach the coherent evolution epoch after matter-radiation 
equality, they experience a scale-dependent shift from the moment 
they re-enter the horizon to the equality epoch. Gravitational in- 
stability is governed by the interplay between radiative pressure 
resistance and gravitational infall. The different duration of modes 
during this period results in a shape dependence on the power spec- 
trum. This shape dependence is determined by the ratio between 
matter and radiation energy densities and sets the location of th e 
matter-radiation equality in the time coordinate jSong et al.ll2010t) . 

One way to estimate cr,, will be to use fitted P00 for each fit- 
ting step. Our measurement is, however, limited at scale of A: < 
0. 1 /; Mpc"' and the contribution to cr,, from ^00 atk >0.1h Mpc"' 
is small but not negligible (~ 10%). Therefore, we calculate cr,, us- 
ing the linear shape of P00 with an amplitude which is estimated at 
each fitting step as follows. 

For each P00 we want to test, we calculate the amplitude fac- 
tor g@(ki, z) defined by 

P@@(k„ z) = gl(k„ z)P@@(ki, ZkJ, (9) 
and constrain the amplitude by calculating a weighted average of 

2;:t:„, {8s(ki,z)/<Ti^(k„z)) 



i:;ZJ/<(k„z) 



Here cr^gikj, z) is given by 



<Tgeik„z) = gl\k„z) 



o-lPl%(ki,z)] 
Pl%(k.,z) 



(10) 



(11) 



and o-[P^f^{ki,Zj)] is given by theoretical estimation in Eq. [5]and 
superscript 'fid' denotes the fiducial model for Fisher matrix analy- 
sis. We would not expect that fractional error of /'©©(fc,, z) is much 
dependent on different fiducial models. The value of cr,, at the best 
fitted power spectra is 2.8/;"' Mpc (the linear theory prediction is 
3.2/7-' Mpc). 
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Pgg(kj) determines the overall amplitude of Pfi,(fc,,/jy), and 
Pggiikj) determines the running of Pfi,(A:,,/jy) in the fi direction. 
These distinct contribution allows us to separate information of 
Pggikj) and Pse(kj) from 5 different fi bins at each kj bin. We find 
these Pgg(ki) and Pe@{ki) by minimizing 



imas 5 5 



X Co\ i,\{ki)[P^^{ki,n^) - /'fl,(fc„yu,)] , 



(12) 



where fc,-„„, = 0.03 /iMpc"' and yt,,„„ = 0.09 /iMpc"'. Off diagonal 
elements of the covariance matrix are nearly negligible and those 
diagonal elements are written as 



Cov-i(A:,) = — . — 1 

a[P,^(k„ti„)f 



(13) 



We present the difference between P(n{ki,iip) (Kaiser effect) and 
Pf,t(kj,iip) (including dispersion effects) in Fig. [2] With the fiducial 
cr,h = 0.24, only one bin of mode k^ = 0.09 /iMpc"' at /j,, = 0.9 
is removed from fitting. Altough this fitting procedure leads to cor- 
relations among different k bins through cr,., those are minimally 
correlated and the results shown Fig. [3] are consistent with theoret- 
ical predictions. 



3 RESULTS AND DISCUSSION 

Velocity-velocity spectra are remarkably well extracted out of 
measured spectra in redshift space at scales k = 0.03, 0.05, 
and 0.07 /iMpc"', and relatively well extracted at scale k = 
0.09 /iMpc"' with more conservative confidence on the threshold 
limit. Filled blue contours in Fig.[3]represent fitted value of Pggiki) 
and Pggikj), and unfilled black contours represent estimation from 
theory with central values given by simulation. For scales from 
k = 0.03 to 0.07 /jMpc"', the decomposed Peeiki) though our fit- 
ting strategy is trustable, which suggests that the few assumptions 
made in this paper are valid for those scales: 

• The assumption of perfect correlation between halo distribu- 
tion and velocity field is correct. The agreement of (A:, ) between 
fitted and true values supports our assumption of r;, ~ 1 indirectly. 

• Dispersion effect is reasonably modelled at scales within our 
confidence limits, which enables us to extract Pggikj) in model in- 
dependent way using estimated cr,,. 

For k = 0.09 /jMpc"', more conservative threshold limits 
should be applied to remove non-linear supression. In Fig. [4] we 
present best fit P@@(ki) with different threshold limits of cTth = 0.24 
(left panel) and cr th =0.18 (right panel). With iXth = 0.24, only one 
bin at /jj = 0.9 is removed. Shown in Fig. |2l extra suppression is 
also observed at fjj = 0.7 bin atk = 0.09 A Mpc"' which can be re- 
moved by more conservative bound cth, =0.18. Shown in the right 
panel of Fig. |4l true ^00 (A:, ) is restored at the cost of weakened 
constraint. 

Theoretical estimation from Fisher matrix analysis is an op- 
timistic bound on errors. It is noticeable that measured varinaces 
(filled blue contours in Fig. [3]( are consistent with estimated vari- 
ances (unfilUed black contours in Fig.[3]l, which assures us that our 
method is optimized extraction of Pggikj) for the given simulation 
specification. 
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Figure 4. Contour plots are shown for decomposed Pgg{kj) and P^^ikj) 
at k = 0.09 /iMpc"' with (T,h = 0.24 (left panel) and 0.18 (right panel). 
Unfilled black and filled blue contours represent the same in Fig.|3] 



4 CONCLUSION 

We propose a statistical tool to decompose Pgg(k) and Peeik) prac- 
tically out of redshift distortion maps, with a few assumptions: 

1) perfect correlation between density and velocity fluctuations, 

2) confidence on theoretical prediction of velocity dispersion ef- 
fect within threshold limit. The results show that the true value of 
velocity-velocity spectra up to = 0.07 /iMpc"' are successfully 
recovered using theoretical dispersion effect. The detectability is 
extendable up to ~ 0.1 /iMpc"' with more conservative threshold 
limit at the cost of weakened constraint. We find that the theoreti- 
cal dispersion effect can be estimated from P@@ik) parameters us- 
ing weighted average at A; < 0.1 hMpc^K In linear regime, Peeik) 
is well-measured with this estimated cr,, as much as with the true 
fixed cr,- of the simulation. 

We find that the biased measurement of Pee(k) is mainly 
caused by the unpredictable non-linear supression effect at k > 
O.lh Mpc"' . The detect ability limit in sca le can be extended by pa- 
rameterizing this effect jTang et al]|2010l) . but we scope our range 
of interest in linear regime in this paper. 
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